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Abstract 

In the context of the minimal supersymmetric standard model, nonzero neu- 
trino masses and mixing can be generated through renormalizable lepton num- 
ber (and thus R-parity) violating operators. It is examined whether neutrino 
mass matrices from tree and one-loop contributions can account for two mass- 
squared differences and mixing angles that explain current experimental data. 
By accommodating, in particular, the solar and atmospheric neutrino data, 
we find interesting restrictions not only on the free parameters of the theory, 
such as lepton number violating couplings and soft-parameters, but also on 
the oscillation parameters of atmospheric neutrinos. 

PACS number (s): 12.60. Jv, 14.60.Pq 
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I. INTRODUCTION 



There exists some evidence for nonzero neutrino masses and mixing. The observations 
of the solar neutrino deficit have been indicating neutrino oscillation . The resonant 
conversion of v e inside the Sun [[| would provide the most favorable explanation for the 
solar neutrino data. Several neutrino experiments have also observed deficit in the atmo- 
spheric neutrino flux [|]]. The evidence for atmospheric neutrino oscillations was recently 
presented by the Super-Kamiokande group |J, which favors v il -v T (or v s ) oscillation. The 
v^-Ve oscillation interpretation can be ruled out by CHOOZ experiment 0. The last one is 
laboratory evidence for the neutrino oscillation — > v e coming from the LSND experiment 
0. This evidence is recently being challenged by KARMEN experiment ||, and is to be 
checked in the near future. The above neutrino data are known to require three distinct 
mass-squared differences and mixing angles; 



Am s 2 ol ~ (4 - 10) x KT 6 eV 2 
sin 2 2# so i ~ (0.12 - 1.2) x 10~ 2 

Am 2 tm ~ (0.5 - 6) x 10" 3 eV 2 
sin 2 2fl atm ~ (0.82-1) 

Am LSND - (0-3 - 2.2) eV 2 
sin 2 2^ LSN D ^ (0.1-4) x 10- 2 




(1) 
(2) 
(3) 



With ordinary three neutrinos, any two of the mass-squared differences in the above equa- 
tions can be obtained: that is, those corresponding to (i) solar and atmospheric (S+A), (ii) 
solar and LSND (S+L), or (iii) atmospheric and LSND (A+L) neutrino data. For (S+A), the 
LSND result has to be disregarded. In the case of (S+L) or (A+L), the presence of a sterile 
neutrino is necessary for the explanation of the atmospheric or solar neutrino experiment, 
respectively. 

One of the desirable features of the supersymmetric extension of the standard model 
would be the generation of small neutrino masses within its context, as the supersymmetric 
standard model with the minimal particle content (MSSM) allows for the lepton (L) and 
baryon number (B) violating operators. In order to ensure the longevity of a proton, one 
usually assumes the conservation of R-parity, forbidding both (renormalizable) B and L 
violating operators. As a consequence, the lightest supersymmetric particle (LSP) is stable 
and thus cold dark matter of the universe may consist of neutral LSP's. However, there is 
no obvious theoretical reason why R-parity needs to be conserved, or why both B and L 
conservation have to be imposed. L-violation would be present in the MSSM and it may be 
the origin of nonzero neutrino masses and mixing that explain current experimental data, 
while proton stability is ensured by B conservation alone. 

The L-violating operators in the MSSM are 

HiLiH 2 , X' ijk LiQjDl , and X^LiLjEl . (4) 
As is well-known [0]], ordinary neutrinos can obtain nonzero masses in tree-level via nonzero 



vacuum expectation values (VEVs) of sneutrinos, as well as at the one-loop level through 
squark or slepton exchanges. Recently, there have been many works studying neutrino 
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phenomenology in the context of R-parity breaking supersymmetric models [ 12] . Typically, 
the tree mass is much larger than the loop mass. The key observation we wish to emphasize 
is that it is, however, possible to find the soft supersymmetry breaking parameter space for 
which the tree mass is rather close to the loop mass, and the solutions to the actual neutrino 
problems can be provided. This was first recognized by Hempfling in Ref. [12|] where a scatter 
plot study of the supersymmetric grand unification model allowing only bilinear operators 
shows that the solar and atmospheric neutrino data can be accounted for. 

In this paper, assuming the presence of trilinear L-violating terms, we will examine how 
the soft parameter space viable for neutrino physics is constrained depending on the choice 
of the L-violating couplings and tan/5. Furthermore, we will find that there are certain 
correlated patterns among soft parameters and predicted masses and mixing for given L- 
violating couplings and tan f3. 



II. PATTERNS OF NEUTRINO MASS MATRICES 

Before investigating neutrino masses arising in the L-violating MSSM, let us first discuss 
the patterns of mass matrices required for each case (i), (ii) or (iii). 

(i) (S+A) 

(a) hierarchical neutrino structure: The R-parity violation may generate at least two nonzero 



m„ 3 , and thus Am^ 



m 2 U2 and Amf tm 



m 



mass eigenvalues satisfying < m U2 < 
required for the explanation of the solar and atmospheric neutrino data. Then, we need a 
mass matrix in the v T directions which yields a large mixing, sin 2 2$^ ~ 1, and a rather 
small ratio between two eigenvalues, 



X = m u Jm U2 = (7-40) 



(5) 



as can be obtained from Eqs.(l) and (2). In this case, the components along the v e direction 
should be able to explain the mixing of v e with v^ r reproducing sin 2 26^ given in Eq. (1). 
(b) degenerate neutrino structure: Another way to accommodate the solar and atmospheric 
neutrino data is to have almost degenerate three neutrinos. This could be achieved in our 
scheme if the amount of L- violation along the e, \x and r directions are the same. Then, 
the solar neutrino data could be explained by a large mixing resonant conversion or vacuum 
oscillation effect [[J which will be disregarded in this paper. 

(ii) (S+L) 



For the explanation of the atmospheric neutrino oscillation in addition to the solar neutrino 



and LSND data, a sterile neutrino almost maximally mixed with has to be invoked |L4 

u s ) basis, 



Let us denote the 4x4 neutrino mass matrix by in 



( rn e 



\ 



^epu TTtfifi TYlfir TTlfis 
Tfl er ?Tl [jt Tfi TT Tfl TS 

\ ?T^es TJlfis TTlss 



(6) 
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where the components mj S come from a certain origin beyond the MSSM. A natural way to 
obtain a large mixing required by the atmospheric neutrino experiment is to have the almost 
Dirac structure with m^ s 3> m es ,m TS . There exist two possibilities to realize this (S+L). 

(a) m Ve1 m VT 3> m v ,m Vg : This is the case where the largest mass scale for Am 2 SND is 
determined by the (e, r) block of rriij which is larger than the (/x, s) block. As the ^ e -z/ M 
oscillation explains the LSND experiment, the solar neutrino data is explained by the v e -v T 
oscillation. Therefore, it is required that v e and v T are almost degenerate: m Ve ~ m Vr ~ 1 eV, 
and m Ur —rn Ue ~ 10~ 6 (10 -10 ) for the MSW (vacuum) solution ||. This extreme degeneracy 
is, however, hard to achieve in our scheme as will become clear in Section IV. 

(b) m u ^,m Va ^> m Ue ,m UT : This possibility is to have v er lighter than u^, for which the 
solar neutrino data can be explained by small mixing in Eq. (1). Then, the mass-squared 
differences for the atmospheric and LSND data are determined by 

Am LSND- m ls, Am atm-H S ( m W+ m »)' ( 7 ) 

The mass-squared difference Am 2 tm for the v^— v a oscillation can be different from that for 
the v }l —v T oscillation due to matter effects in the Earth. Recent analysis ]13|] based on the 



Super-Kamiokande data shows that Am 2 tm = (0.2 — 1) x 10 2 eV 2 for the v^-v s oscillation 
which is a bit shifted up compared with the value in Eq. (0). For this value of Am 2 tm and 
A^lsnd in Eq. (3), we get 

~ (0.55 - 1.5) eV 
m^ + m ss ~ (1.3- 18) x KT 3 eV . (8) 

Since the mixing angle required for the explanation of the LSND result is given by #lsnd — 
rries/m^s, one needs m es ~ (2.4 — 5.5) x 10~ 2 eV. 

Let us now estimate the sizes of m^- for active neutrinos coming from R-parity violation. 
If the lepton number violation in the /i direction is suppressed so that m^,m^ e ,m^ T <C 
m ss , one needs to generate only the v e —v r oscillation for the solar neutrino for suitable 
values of m eT and m TT which can be rather trivially obtained in terms of the lepton number 
violations in the direction of e and r. This case requires an explanation for the origin of 
the mass scale m ss ~ 10 -2 eV. It would be more natural to have m ss w as in Ref. ITU 



Vgs ' v iu \j v . it vvuuiu vj\^ i±iv7i v natui ai iu nave 'it's. 

and thus m MA1 ^> m ss . In this case, one needs m M(U — (1-3 — 18) x 10~ 3 eV. Under the 



assumption that m TS /m^ s <C 1, the solar neutrino data can be explained if m TT = yAm 2 ol 
and m eT /m TT — m^ T m es /m TT m^ s = 9 m \. From this, the component m^ T has to be restricted 



as m^ T < J Am 2 ol ^ so i/^LSND- Therefore, as in the previous case (i), the L- violation in the 
direction of /i and r has to generate not so large ratios: 

m w /m TT ~ (0.41 - 9) , m^ T /m TT < 3.5 . (9) 

For this calculation, we have used the value of m^ s varying from 0.55 eV to 1.5 eV and the 
largest value of m^ T . 

(iii) (A+L) 
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This is the most popular case with four neutrinos [HJ. The solar neutrino problem is 



solved by v e —u a mixing, and the atmospheric neutrino problem by maximally mixed ^-/x T 
oscillation. Typically, one would expect that u e jS are lighter than z/ MT . Combined with 
the LSND result, it is then required that the two mass eigenvalues satisfy m V2 ~ m Vz ~ 
(0.55 — 1.5) eV (therefore \ ~ 1) and m 2 ^ — m 2 2 = Am^ tm . The LSND mixing angle requires 
also m^/m^ or m eT /m^ T ~ (2.4 — 5.5) x 10~ 2 . In the scheme of generating the neutrino 
masses from the R-parity violating couplings, it is again difficult to produce a small number: 
(m U3 — m U2 )/(m U3 + m V2 ) = Am^ tm /4m^ ~ 10 -3 . Recall that it is also possible to have 
almost degenerate z/ e s with mass around 1 eV, and lighter jT . In this case, one would need 
to fine-tune some parameters in the v e -v s mass matrix to achieve small mixing and a very 
small mass difference, ~ or 10~ 10 eV. 



III. MSSM WITH R-PARITY AND L VIOLATION 



Our framework is the conventional MSSM [16| in which soft supersymmetry breaking 



parameters arise from the gravitational coupling to the hidden sector and are assumed to 
be universal at the grand unification (GUT) scale where the three gauge couplings meet. In 
this framework, there are five independent parameters; the scalar mass mo, the trilinear soft 
parameter A , the gaugino mass m^, as well as tan/? and the sign of /i. Here fi represents 
the mass parameter of the bilinear operator H1H2. 

Without loss of generality, one can take only dimensionless L-violating couplings A', A 
at the GUT scale. Recall that, in models where all the L-violating terms appear generically, 
the universality condition allows us to redefine the superfields Li and Hi in such a way that 
//j = Dj = at the GUT scale. Here Di is a dimension-two soft parameter corresponding to 
the bilinear operator L1H2. Even if fii and Di are zero at the GUT scale, upon renormalizing 
the soft SUSY breaking parameters and the L-violating couplings down to the weak scale, 
the universality condition breaks down and nonzero /ij, Di and vn 2 LiHi are generated. For 
the calculation of the renormalization group equations (RGEs), it is convenient to use the 



basis [17 where the /ij terms are continuously rotated away in the following approximate 



manner: 



Fx + —Li — > H x (10) 
Li H\ — > Li . 

This definition of new basis is valid up to the leading order of A' and A, which can be 
justified if they are small enough. The merit of this basis is that the RGEs of A' and A do 
not mix each other as shown in Appendix. Nonvanishing Di and rn? L . Hi at the weak scale 
induce nonzero VEVs of sneutrinos |jTT|. Keeping only the leading terms of the L-violating 
soft-parameters, one finds the sneutrino VEVs [|l7j]: 

(p.) A^ + mWi ( s 

m 2 Li + Mf cos 2/3/2' K ' 

The tree-level neutrino mass matrix is then given by 
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(mT) = 1 {Dl){&3) 5 ^- (12) 

v '»j mi/2/47raG[/T + 1> 2 sinp//i 

Given the L-violating couplings A' and A at the weak scale, the neutrinos acquire radiative 
masses whose matrix elements are given by 



m 



ioo P \ ^ 3 A dk m dl 
'ij IQtc 2 fhj 
1 A 



+ i6^^r L [^m + (^j)\, (13) 

where m^e denote the down-type quark and charged lepton masses, respectively, Ad >e denote 
the mixing masses of the corresponding squarks or sleptons, and m 2 are the functions of the 
squark or slepton mass eigenvalues m lj2 given by m 2 = (m 2 — m 2 )/ In (m^/m 2 ,) for each 
generation. Here, k, I are generation indices. The full neutrino mass matrix is now given 
by (m„) = (m}J ee ) + (m}° op ), which will be computed with specifying the input values: the 
L-violating couplings A', A, the soft parameters mo, A , m\/2 (fixed at the GUT scale), as well 
as tan f3 and the sign of \i. For the computation, we took the ranges of the soft parameters 
as follows: mo from 100 GeV to 1 TeV, the absolute value of A from 100 GeV to 3mo, and 
the absolute value of mm from 100 GeV to 1 TeV, fixing the sign of \i to be positive. Note 
that changing the signs of /i, A and m\j% simultaneously yields the same results. We also 
fixed the top quark mass to be 175 GeV, and the strong coupling constant a s [Mz) to be 
0.118. 

Let us now specify more on the required L- violation. Since there are too many L-violating 
couplings, it would be impossible to draw sensible conclusions allowing such arbitrariness in 
explaining the neutrino oscillation parameters in Eqs. (l)-(3). Our basic assumption in this 
regard is to take the couplings with a "natural" hierarchy: the L-violating Yukawa couplings 
have a hierarchical structure similar to the corresponding quark or lepton mass matrix. In 
other words, e.g., A^- fe oc h^ k where hj k is the Yukawa matrix for the down- type quarks before 
diagonalization. This will be a consequence of models that explain the fermion hierarchies 
in terms of flavor symmetry [ 18[| . In such a scheme, one expects that A^ 33 and Aj33 are the 



largest components and thus give the leading contribution to the neutrino masses. The 
subleading contribution comes from the couplings, e.g., A- 32 or A- 23 . Their contribution to 
the component {m v )ij = rriij would be suppressed by the factor of (A^ 32 Aj- 23 /A^ 33 Aj- 33 )(m s /mb) 
compared with the leading contribution. It is to be understood that this expression is taken 
after diagonalization of the quark mass matrices. It was realized recently by Drees et.al. 
|19| that the ratio m s /mb could account for the ratio x" 1 if -^32^23 ~ ^33^-33- But, for 
the couplings with the natural hierarchy, it is expected that A^ 32 Aj 23 /A^ 33 A^ 33 < m s /m b , and 
thus the contribution of the smaller couplings is suppressed by a factor (m s /mb) 2 . This 
is too small to account for the size of various components of the neutrino mass matrices 
leading to the case (i), (ii) or (iii). Therefore, the couplings other than the (33)- components, 
K33 ( or \2>z)> can be neglected in our discussion. Now, the question one may ask first 
is whether the phenomenologically desirable neutrino mass matrices can be obtained by 
taking a minimal set of the L-violating couplings, that is, the trilinear couplings other than 
A^ = A^ 33 are negligibly small. After answering this, we will examine the effect of the 
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couplings Aj33 = Xi(i 3). We concentrate on the couplings of X' 2 3 or A2 since the coupling 
\[ or Ai (smaller than others) can be adjusted to reproduce the desirable mixing of v e to 
other neutrinos as already discussed. 



IV. NEUTRINO MASSES FROM R-PARITY VIOLATION 

Our task is now to examine whether the neutrino masses coming from the R-parity 
violation in the MSSM can reproduce the mass matrix patterns studied in Section II. Before 
presenting our main results, it is useful to look into some qualitative features of the neutrino 
masses coming from R-parity violation. First thing to note is that typically the tree mass 
is much larger than the loop mass. This can be seen by taking a crude approximation 
to integrate the RGEs in the Appendix. That is, taking the constant coefficients in the 
RGEs for D{ and rn? L Hl , one can easily obtain their values at the weak scale which yield the 
sneutrino VEVs 

{n) « vt Mz ai (X% + hXiK) (14) 

where ai, bi are the parameters encoding the RGE effects for the soft masses and couplings. 
Therefore, the tree and loop mass can be written as 

(m* ree ) . . « KoaittjiX'ihb + hXihr) (X'jh + bjXjh T ) 

(m'r)^ « ^i(K^hl + bijXiXjhl) (15) 

where k ~ (31n(M GC / T /Mz)/87r 2 ) 2 (Mf cos 2 /3/mi/ 2 ) and Ki ~ 3u 2 /87r 2 m. Generically, one 
would find m trec / m lo °v> ^ ko/ki > 10 2 for m/m^ ^ 0.1. The second observation is that, 
for A' ^> Aj, the tree mass can be aligned with the loop mass; that is, a« ~ aj for i ^ j. 
The alignment can be weakened as tan/5 and Aj become larger. For a large tan/3, the 
misalignment of, e.g., A' 3 /X' 3 — A' 2 /X' 2 can be amplified by the RGE effect (and thus a 3 — a 2 
also) since 

J / AI M \ 

a h T A T . (16) 

For large Aj, it is obvious to have a large misalignment since A3 = due to the antisymmetry 
of Xijk- 

The above properties play important roles in the case (S+A), for which one needs two 
distinctive nonzero mass eigenvalues. Since the tree mass is much larger than the loop mass, 
the loop mass can fit better the solar neutrino mass scale with vacuum oscillation (requiring 
Am 2 ol ~ 10 _10 eV 2 ) while the tree mass gives rise to the atmospheric neutrino mass scale. 
This is basically why one finds more soft parameter space for the vacuum oscillation solution 
to the solar neutrino problem than for the MSW solution in the scatter plot study by 




Hempfling []T5|. The vacuum oscillation of solar neutrinos requires a large mixing between 
v e and v T in our scheme since the large mixing for vacuum oscillation implies X[ ~ A23. 
This might be in conflict with the Superkamiokande data as well as the CHOOZ result 
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|5|||. Therefore, we prefer the small mixing MSW solution to the solar neutrino problem as 
mentioned before. It is now clear that, in order to get the neutrino mass hierarchy viable for 
the atmospheric neutrino oscillation and the MSW solar neutrino conversion, the tree mass 
has to be suppressed requiring a, <C 1. This can occur for certain soft parameters which 
admit a cancellation in (z/j), that is, rn 2 LiHl /Di ~ — tan/3. Then, it could be that the ratio 
of the tree and loop mass is responsible for the two distinctive neutrino masses differing by 
the factor x = (7 - 40). 

With the minimal number of L- violating couplings (or, Aj <C A^), the mere suppression 
of the tree mass is not enough because the tree and loop masses are almost aligned as 
discussed above (that is, (z>j) oc A^), which renders the second eigenvalue m U2 close to zero. 
However, it turns out that a partial cancellation between the tree and loop mass can result 
in the breakdown of the alignment and the production of the desirable value of \. In other 
words, it is required that |m* ree | ~ |m|, oop | and they have opposite sign. In order to quantify 
these statements, let us give some examples and calculate the amount of the misalignment 

measured by f = \\' 2 (v 3 ) ~ Ki^} |/|A'| | (v) | where |A'| = y ' \' 2 2 + A' 3 2 , etc. For tan/3 = 5, 
m = 201 GeV, m\/2 = —109 GeV, and A = —550 GeV, the muon and tau neutrino mass 
matrix is found to be 

m " /6V U-68 0.69j + U.65 0.65 V 1 [ ' 

where the first (second) one is tree (loop) contribution. As is obvious from Eq. fll7|) , the tree 
and loop masses are aligned very closely (£ = 7 x 10 -3 ) and the cancellation of the largest 
digit results in x ~ 38 and sin 2 29 « 0.9. For tan/3 = 40, taking m = 20 GeV, mi/ 2 = —110 
GeV, Aq = —250 GeV, one finds the mass matrix with x ~ 16 and sin 2 29 ~ 0.9; 

/0.69 1.22 \ /4.01 3.87\ f . 

m »/ eV *-[l.22 2.14J + U.87 3.73J (18) 

which shows that the alignment is badly broken (£ = 0.1). Still, one needs a cancellation 
to get the right value for the mass ratio, in particular. It turns out that negative values 
of mi/2 and A are necessary to yield opposite signs for the tree and loop masses. Now 
that the alignment becomes severer for smaller tan f3, the soft parameters have to be tuned 
more precisely for smaller tan/3. In the numerical computation, we calculated the ratio x 
of two mass eigenvalues and the atmospheric neutrino mixing sin 2 29 for some sample set 
of tan/3 and the ratios between the couplings A' 23 , scanning the soft parameters up to the 
smallest digit above point. The results of the computations are presented in TABLE I and 
II. Notice that the neutrino mass matrix is proportional to the overall scale of A'23, an d 
thus the ratio x an d the mixing sin 2 29 depend only on the ratio of two couplings A 3 /A 2 . As 
expected, the number of the desirable soft parameters gets smaller for smaller tan/3, and 
we find no parameter space for tan (3 < 4. Furthermore, the values of x an d sin 2 29 depend 
very sensitively on the soft parameters, and the acceptable soft parameters are scattered for 
a small tan f3. For a large tan /3, x becomes a slowly varying function of the soft parameters, 
and one can isolate some finite region for given values of x an d sin 2 29. 

Unexpectedly, some patterns have emerged for values of the ratio x an d the mixing 
angle that are realized in our scheme depending on tan (3 and the L- violating couplings, as 
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presented in TABLE I. For illustration, we took some variations of A3 larger than A' 2 . Similar 
results can be obtained also in the cases with A' 2 larger than A 3 . For given tan/3 and A 23 , x 
and sin 2 29 are found to lie inside some restricted ranges, and to be correlated in the way that 
larger mass ratio corresponds to larger mixing for A 2 = \' 3 , or to smaller mixing for A 2 < A 3 . 
The correlation becomes weaker for larger tan j3. From TABLE I, one sees that the maximal 
mixing is easily realized with A 3 /A 2 ~ 2. We found practically no desirable parameter space 
for \y X' 2 > 5 which appears inconsistent with large sin 2 29. The eigenvalue m U2 depends 
on mo in the way that larger tuq produces smaller neutrino mass as implied by the formulae 
(|TT|)-(|T^). The right values of the L- violating couplings can be obtained by the rescaling: 
A 23 = 10~ 4 (Am 2 ol /m 2 2 ) where m V2 is a value given in TABLE I. For instances, the actual 
values one needs are A 2 = A' 3 = (1(T 3 - 1(T 4 ) for tan/3 = 5, A' 2 = A 3 = (10~ 4 - KT 5 ) for 
tan (3 = 20 and A 2 = A' 3 = (3 x 10~ 5 — 3 x 10~ 6 ) for tan j3 = 40. As a consequence of partial 
cancellation between the tree and loop mass, the size of the L-violating couplings yielding 
the phenomenologically viable neutrino masses becomes slightly larger than the values for 
which the loop mass yields the atmospheric neutrino mass scale, say 5 x 10~ 2 eV. 

It is also amusing to find that there are preferable ranges of the acceptable soft parame- 
ters: larger \A \ and smaller |mi/ 2 | are preferred. To be specific, we show in TABLE II the 
ranges of soft parameters within which viable neutrino masses and mixing can be realized 
with A' 2 = A 3 = 10~ 4 . For a given m , the values of A and mi/ 2 should reside in between 
two end values shown in the table. There is a correlation between Aq and m^: larger \Aq\ 
corresponds to larger |m!/ 2 |. A similar pattern follows even with a slight variation of two 
couplings A' 2 3 such as given in TABLE I. 

If the coupling A2 is comparable to A 23 , then the above mentioned properties are signif- 
icantly modified as anticipated in the beginning of this section. We find that x an d sin 2 29 
become very slowly varying functions of the soft parameters even for a small tan/3, as a 
consequence of a large misalignment between A- and (pi). But one still needs m* ree ~ m° op 
requiring a cancellation in the sneutrino VEVs, which essentially restricts the soft parameter 
space. To show this explicitly, let us take two examples for tan/3 = 5. Taking m = 500 
GeV, mi/2 = —100 GeV and Aq = —780 GeV, one finds the mass matrix in which the tree 
mass is larger than the loop mass as in the above minimal cases: 



giving x ~ 38 and sin 2 29 « 0.89. There is a small contribution (0.003 eV) to the m l ° ° p 
from the A 2 coupling. The other one is of different type with dominant loop contribution: 
for m = 500 GeV, m 1/2 = 400 GeV and A = 100 GeV, one gets 



which gives x ~ 20 and sin 2 29 ~ 0.95. In this case, the A2 contribution (0.009 eV) is more 
significant than before. In both examples, any of the tree and loop contributions cannot be 
neglected in order to produce the right values of x an d sin 2 29. In TABLE III, we illustrate 
some set of soft parameter ranges yielding the right values of the ratio x an d the mixing 
angle, taking A' 2 = A' 3 = A2. The ranges of A and rai/2 for given mo are shown. In TABLE 




(19) 




(20) 
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IV, we take some variations of A2/A23. Compared to the previous case, one finds more 
parameter space open up for larger A2. To have the right value of m U2 , we need a bit smaller 
values for the couplings A' than before: that is, A' 23 = A 2 — (1CT 4 — 3 x 10~ 5 ) for tan/? = 5 
and A 23 = A 2 ~ (10~ 5 - 2 x 10~ 6 ) for tan/3 = 40'. From TABLE III and IV, one can also 
see that larger tan/3 are still needed to destroy the alignment in a sufficient amount even 
with sizable A2. It can be found that a large soft parameter space opens up for tan/3 < 2 
only when A2 ^ 4A 2 . Contrary to the case with negligible A2, almost all ranges of x> sin 2 29 
can be realized as shown in TABLE III. 

Now let us turn to the other cases. For the case (S+L), it is important to realize 
that one needs not to generate two nonzero eigenvalues in the neutrino mass matrix along 
the Up, v T directions. Furthermore, the ratios m^ T /m rr and m^/m rr (|9|) are determined 
roughly by the input values A 2 /A' 3 and (A 2 /A 3 ) 2 , respectively. Therefore, it is required from 
the previous discussion, A' 2 /A 3 ~ (0.64 — 3), and the value of A' 3 to be determined by the 

condition m v TT = J Am 2 ol ~ 10 _3 eV. Generically, the tree mass is much larger than the loop 
mass. Taking A^ 33 = 10 -4 , the tree mass can be as large as 100 eV for a small tan/3 ||12|| . 
Therefore, for a small tan/3, one needs A 3 ^ 10~ 7 . This value can be as large as 10 -4 when 
the tree mass is suppressed as we discuss above. Since Aj cannot give rise to the nonzero 
component vn v TT , taking the largest components being A' 2 ~ A' 3 as above is the best way to 
explain m v m ~ m^ T ~ m u TT . 

Let us finally comment on the case (A+L). In order to achieve (m U3 —m V2 )/ (m V3 +m U2 ) ~ 
10~ 3 , much finer tuning of the soft parameters is required and thus it is very difficult to be 
realized in our scheme. 

V. CONCLUSIONS 

In conclusion, we have examined the possibility of obtaining the realistic neutrino masses 
and mixing in the context of the R-parity violating minimal supersymmetric standard model. 
We have assumed an ultraviolet theory which has generic L-violating Yukawa couplings in 
the basis where the L-violating bilinear terms are rotated away. The L-violating Yukawa 
couplings induce the L-violating bilinear soft terms through the renormalization group evo- 
lution which takes a simple form in the basis (valid for small L-violating couplings) as shown 
in Appendix. Analyzing the neutrino masses arising both from the sneutrino vacuum expec- 
tation values generated by the L-violating soft terms and from the loop corrections through 
squark or slepton exchange, we found restrictions on the soft parameters, the L-violating 
couplings and tan /3, under which realistic neutrino mass matrices can be obtained. 

With three ordinary neutrinos, one can account for any two of the three distinct mass- 
squared differences required by the solar, atmospheric and LSND neutrino data. First, we 
have discussed the phenomenological mass matrices along the z/ M , v T directions which are 
required by the data. If the solar neutrino and LSND data are due to the active neutrinos, 
and a sterile neutrino is introduced to explain the atmospheric data, then it is not necessary 
to have two distinct mass eigenvalues for the mass matrix in the u^, v T directions but its 
components should not differ by a factor of more than a few. This can be easily achieved 
if A^ 33 are the largest L-violating couplings and A' 233 ~ A' 333 . In the case of solving the 
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atmospheric neutrino and LSND data, one needs to generate two almost degenerate mass 
eigenvalues for v T with a very small mass difference. This case can hardly be realizable 
in our scheme. 

Most nontrivial and realistic case is to accommodate the solar and atmospheric data 
within the ordinary context of three active neutrinos (disregarding the LSND data). In this 
case, one needs to generate two distinct mass eigenvalues for v T whose ratio, Xi should 
reside roughly between 7 and 40, and the mixing angle sin 2 20 atm ~ 0.82. 

Under the assumption of the natural Yukawa hierarchy in the L-violating couplings, we 
have argued that the relevant contribution to the phenomenological neutrino mass matrices 
comes from the components A^ 33 and A^. With the minimal choice of the L-violating 
couplings (namely, other than A^ 33 are negligible), we needed not only the suppression of 
tree mass (that is, tan/3 ~ — m 2 L . Hl /Di), but also some partial cancellation between the tree 
and loop mass. This basically strongly constrains the soft parameter space. Our study have 
shown that the realistic neutrino masses and mixing prefer a large trilinear soft parameter A 
and a small gaugino mass m^. The desirable soft parameter space becomes more restricted 
and finer tuned for smaller tan/5, and thus a reasonable parameter space can be found 
only for fairly large tan/3, say, tan/3 > 40. We have found indeed no parameter space for 
tan (3 ^ 4. The large mixing explaining the atmospheric neutrino data requires that A 233 and 
A 333 should not differ by more than a factor of 5. Interestingly, the mass ratio \ = m u 3 / m u 2 
and the atmospheric neutrino mixing sin 2 29 are found to be restricted in a certain range 
and correlated in the way that smaller ratio \ has smaller or larger mixing depending on the 
values of A 2 3 . This tendency becomes also weaker for larger tan j3 and larger A233. It appears 
more reasonable to have a sizable A233, for which the experimental quantities of neutrino 
oscillations become very stable under the variation of the soft parameters, and thus there 
exist fairly extended regions of parameters fitting into the experimental data. Still, there 
exist significant constraints on the soft parameters coming from the fact that the tree mass 
has to be suppressed as above. More soft parameter space can be found for larger tan/3 and 

A233- 



APPENDIX 



Renormalization group equations for the lepton number violating parameters in the basis 
where LiH\ terms are rotated away in the superpotential. 

16vr 2 ^ = AK« + hi + 3hl - 7 -gl - 3, 2 - fgl) (21) 

16vr 2 ^ = HShl - 3gf - 3^ 2 ) (22) 

d A' 7 Ifi 

16vr 2 ^ = A[(5 lZ hl + h\ + 9h" b - L g l _ 3g l _ J^g) + Ai (2h b h T ) (23) 

7 16 

+ 2A^ 3 /i t At + Mi + 2h b A b + -g 2 ^ + 3# 2 M 2 + y^M 3 ) 

dA 

167r2 ^T = M$h 2 T - 3</ 2 - 3<? 2 ) + A[{Qh b h T ) + A,(6/i r A r + §g\M x + 6^ 2 M 2 ) (24) 
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16tH 



.dm 



dt 



,dD 



m 2 LiHl ^S i3 h 2 T + h 2 T + 3h 2 b ) - 6A' t A b - 2A,A T 
6A-/i b (m|. + m 2 Qs + m 2 D3 ) - 2A i /i T (m|. + m| 3 + m|J 



1671 df = Dli * 1 ~ 91 ~ %2) " ^ 6A ' ih + 2AlK) 



(25) 
(26) 



ACKNOWLEDGMENTS 

EJC would like to thank B. de Carlos for useful communications. UWL thanks KIAS 
for the kind hospitality during his visit, and BSRI of Mokpo National University. EJC is 
supported by Non-Directed Research Fund of Korea Research Foundation, 1996. 



12 



REFERENCES 



[1] J.N. Bahcall and M.H. Pinsonneault, Rev. Mod. Phys. 67, 781 (1995); J. N. Bahcall, S. 

Basu and M. H. Pinsonneault, Phys. Lett. B 433, 1 (1998). 
[2] Homestake Collaboration, B. T. Cleveland et.al., Nucl. Phys. B (Proc. Suppl.) 38, 

47 (1995); Kamiokande Collaboration, Y. Fukuda et.al., Phys. Rev. Lett. 77, 1683 

(1996) ; GALLEX Collaboration, W. Hampel et.al., Phys. Lett. B 388, 384 (1996); 
SAGE Collaboration, J. N. Abdurashitov et.al., Phys. Rev. Lett. 77, 4708 (1996); Super- 
Kamiokande Collaboration, Y. Suzuki, talk at Neutrino 98, Takayama, Japan, June 
1998. 

[3] S.P Mikheyev and A.Yu. Smirnov, Yad. Fiz. 42, 1441 (1985) [Sov. J. Nucl. Phys. 42, 

913 (1985)]; II Nuovo Cimento C 9, 17 (1986); L. Wolfenstein, Phys. Rev. D 17, 2369 

(1978); ibid. 20, 2634 (1979). 
[4] Kamiokande Collaboration, Y. Fukuda et.al., Phys. Lett. B 335, 237 (1994); 1MB 

Collaboration, R. Becker-Szendy et.al., Nucl. Phys. B (Proc. Suppl.) 38, 331 (1995); 

Soudan-2 Collaboration, W. W. M. Allison et.al., Phys. Lett. B 391, 491 (1997). 
[5] Super-Kamiokande Collaboration, Y. Fukuda et.al, Phys. Rev. Lett. 81, 1562 (1998), 

[hcp-cx / 98070031 ; T. Kajita, talk at Neutrino 98, Takayama, Japan, June 1998. 
[6] M. Apollonio et.al, Phys. Lett. B 420, 397 (1998). 

[7] C. Athanassopoulos et.al., Phys. Rev. Lett. 77, 3082 (1996); preprint UCRHEP-E197, 

|nucl-ex / 97090061 ; H. White, talk at Neutrino 98, Takayama, Japan, June 1998. 
[8] B. Zeitnitz, talk at Neutrino 98, Takayama, Japan, June 1998. 

[9] N. Hata and PC Langacker, Phys. Rev. D 56, 6107 (1997); Y. Suzuki, in Ref. @; J. 

N. Bahcall, P. I. Krastev and A. Yu. Smirnov, Phys. Rev. D 58, 096016 (1998). 
[10] M. C. Gonzalez- Garcia, H. Nunokawa, O. Peres, T. Stanev and J. W. F. Valle, preprint 

|hep-ph/97l2~2^ ; V. Barger, T. J. Weiler and K. Whisnant, Phys. Lett. B 427, 97 (1998). 
[11] L. J. Hall and M. Suzuki, Nucl. Phys. B 231, 419 (1984). 

[12] R. Hempfhng, Nucl. Phys. B 478, 3 (1996); B. de Carlos and P. L. White, Phys. Rev. 
D 54, 3427 (1996); A. Yu. Smirnov and F. Vissani, Nucl. Phys. B 460, 37 (1996); H. P. 
Nilles and N. Polonsky, Nucl. Phys. B 499, 33 (1997); F. M. Borzumati, Y. Grossman, 
E. Nardi and Y. Nir, Phys. Lett. B 384, 123 (1996); E. Nardi, Phys. Rev. D 55, 5772 

(1997) ; R. Adhikari and G. Omanovic, preprint |hep-ph/9802"390] 

[13] R. Foot, R. R. Volkas and O. Yasuda, Phys. Rev. D 58, 013006 (1998). 

[14] Q. Y. Liu and A. Yu. Smirnov, Nucl. Phys. B 524, 505 (1998); E. J. Chun, C. W. Kim 
and U. W. Lee, Phys. Rev. D58 093003 (1998), [Eep^ph/9802209[ 

[15] J.T. Peltoniemi and J.W.F. Valle, Nucl. Phys. B 406, 409 (1993); D.O. Caldwell and 
R.N. Mohapatra, Phys. Rev. D 50, 3477 (1994); E. Ma and P. Roy, Phys. Rev. D 52, 
R4780 (1995); E.J. Chun, A. S. Joshipura and A. Yu. Smirnov, Phys. Lett. B 357, 608 
(1995); ibid Phys. Rev. D 54, 4654 (1996); R. Foot and R.R. Volkas, Phys. Rev. D 52, 
6595 (1995); Z. Berezhiani and R.N. Mohapatra, Phys. Rev. D 52, 6607 (1995); E. Ma, 
Mod. Phys. Lett. A 11, 1893 (1996); A.Yu. Smirnov and M. Tanimoto, Phys. Rev. D 
55, 1665 (1997). 

[16] For reviews, H. P. Nilles, Phys. Rept. 110, 1 (1984); R. Arnowitt and P. Nath, preprint 
CTP-TAMU-52/93, NUB-TH-3073-93; M. Drees, preprint |hep-ph/96ll4TJ9l, ; S. P. Mar- 
tin, preprint [hep-ph/9709356 . 



13 



[17] For instance, see Carlos and White in Ref. [12 



[18] E. J. Chun and A. Lukas, Phys. Lett. B 387, 99 (1996); K. Choi, E. J. Chun and H. 

Kim, Phys. Lett. B 394, 89 (1997). 
[19] M. Drees, S. Pakvasa, X. Tata and T. ter Veldhuis, Phys. Rev. D 57, 5335 (1998). 



14 



TABLES 



TABLE I. Ranges of the atmospheric neutrino mixing angle, the ratio of two mass eigenvalues 



corresponding to ^Am^/Am^, and the smallest mass eigenvalue (for solar neutrino) which can 
be realized in the soft parameter space for given tan/3 and L-violating couplings A 23 . A' 2 = 1CT 4 
is taken. 





tan B 


5 




20 


40 




sin 2 26 


82 — 91 

V / • V.- 1 _ v / • • / 1 




82 — 95 


0.82 — 1.0 


A 3 /A' 2 = 1 


m Vz /m V2 


19 - 40 




14 - 40 


7 - 40 




m V2 (eV) 


3 x 10" 5 - 3 x 10~ 3 




10~ 3 - 0.1 


0.03 - 2 




sin 2 26 


0.89 - 1.0 




0.86 - 1.0 


0.82 - 1.0 


^3/ ^2 = ^ 


m u .Jm U2 


40 - 7 




40 - 7 


7 - 40 




m U2 (eV) 


3 x 10" 3 - 10" 2 




0.01 - 0.4 


0.05 - 4 




sin 2 26 


0.82 ~ 0.96 




0.82 - 0.97 


0.82 - 0.99 


^3/^2 = 3 


m U3 /m U2 


16 - 7 




16 - 7 


28-7 




m V2 (eV) 


2 x 10~ 4 ~ 2 x 10~ 2 


5 


x 10~ 3 - 0.8 


0.1 - 5 




sin 2 26 


0.82 - 0.88 




0.82 - 0.89 


0.82 - 0.94 


A' 3 /A' 2 = 4 


m U3 /m U2 


9 - 7 




9 - 7 


9 - 7 




m U2 (eV) 


4 x 10~ 4 - 2 x 10~ 2 




i x 10~ 3 - 1 


0.4 - 2 


TABLE II. 


Illustrated ranges of soft parameters in units of GeV within which the solar and 


atmospheric neutrino data can be 


accommodated by suitable choices of soft parameters in the case 


of A' 2 = A 3 = 10 


4 . The sign of fi is fixed to be positive. 








tan/3 


m 


A 






nT-i/2 


5 


200 


-590 - 


-500 




-150 100 




400 


-1000 - 


-730 




-270 - -100 




600 


-950 - 


-900 




-150 - -100 


20 


200 


-600 - 


-230 




-270 - -100 




500 


-1000 - 


-300 




-430 150 




800 


-1000 - 


-370 




-430 100 


40 


200 


-600 - 


-200 




-300 ~ -100 




500 


-1000 - 


-500 




-500 - -100 




800 


-1000 - 


-240 




-480 100 



TABLE III. Illustrated set of soft parameter ranges yielding realistic values of the mixing and 
the mass ratio in the case of A' 2 = A3 = A2 = 10~ 4 . The sign of /j, is again taken to be positive. 
Soft parameters are in units of GeV and neutrino masses are in units of eV. 



tan/3 


m 


A 




m 1/2 


m U2 


sin 2 26 


m V3 jm V2 




100 


-100 - 


-170 


430 - 570 


(2 - 7) x 10~ 2 


0.82 - 0.96 


7 - 40 


3 


300 


-310 - 


-830 


490 - 210 


(1 ~ 3) x 10~ 2 


0.82 ~ 0.94 


7-40 




500 


-830 - 


-970 


510 - 470 


(0.9 - 1.1) x 10~ 2 


0.82 - 0.88 


25 - 40 


5 


100 
500 


100 - 
-770 - 
100 - 


310 

-990 

530 


160 - 490 

-100 240 

370 - 1000 


(5 - 16) x 10~ 2 
(0.6 - 2) x 10~ 2 
(1 - 3) x 10~ 2 


0.82 - 0.91 
0.82 - 1.0 
0.82 - 1.0 


7 - 40 
7-40 
7-40 
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900 


-100 170 

100 ~ 530 


690 ~ 550 
690 - 1000 


(5 ~ 6) > 
(5 ~ 10) 


< 10 -3 
x 10 -3 


0.99 ~ 1.0 
0.88 ~ 1.0 


37 ~ 40 
14 ~ 40 




300 


-630 910 


-290 430 


2.8 ~ 


3.5 


0.82 ~ 1.0 


7 ~ 12 






con m n 

OoU ~ yiu 


zyu ~ Oil) 


z.o ~ 


6.0 


U.oz ~ l.U 


7 1/1 


40 


500 


-250 330 


-100 130 


0.4 ~ 


0.9 


0.82 ~ 1.0 


7 ~ 19 






930 ~ 1000 


510 ~ 570 


1 ~ 


2 


0.82 ~ 1.0 


7 ~ 9 




900 


-300 590 


-100 230 


0.1 ~ 


0.5 


0.82 ~ 1.0 


7 ~ 40 






100 ~ 280 


110 ~ 210 


0.2 ~ 


0.6 


0.82 ~ 1.0 


7 ~ 40 



TABLE IV. Same as in TABLE III with different set of A2. Here two approximate end values 
of soft parameters A and m 1 / 2 with the corresponding neutrino parameters are shown. 





A 2 = 3 x 10^ 5 


A 2 = 3 x 10~ 4 


tan/3 


m 


A 




A 


^1/2 


3 


100 


-300 120 


120 ~ 240 


100 ~ 160 


380 ~ 1000 




300 


-900 ~ -640 


320 ~ 420 


-100 900 


150 ~ 1000 




900 


none 


none 


-1000 620 


880 ~ 1000 


5 


100 


100 ~ 300 


140 ~ 380 


100 ~ 300 


220 ~ 1000 




500 


-1000 820 


-220 ~ -120 


-1000 100 


-200 ~ 460 






none 


none 


100 ~ 340 


580 ~ 1000 




900 


none 


none 


-1000 ~ -100 


120 ~ 1000 


40 


300 


-300 240 


-120 ~ -100 


-600 320 


-320 ~ -160 






180 ~ 260 


120 ~ 140 


240 ~ 560 


160 ~ 360 




500 


-560 ~ -300 


-240 ~ -120 


-1000 540 


-540 ~ -260 






140 ~ 460 


120 ~ 240 


400 ~ 980 


280 ~ 680 




900 


-1000 420 


-420 ~ -100 


-320 300 


-100 120 






100 ~ 800 


120 ~ 480 


740 ~ 1000 


520 ~ 680 
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